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ABSTRACT 

An initial value approach is used to examine the dynamics of perturbations introduced into 
a vortex under strain. Both the basic vortex considered and the perturbations are taken as fully 
three-dimensional. An explicit solution for the time evolution of the vorticity perturbations is 
given for arbitrary initial vorticity. Analytical solutions for the resulting velocity components 
are found when the initial vorticity is assumed to be localized. For more general initial vorticity 
distributions, the velocity components are determined numerically. It is found that the variation 
in the radial direction of the initial vorticity disturbance is the most important factor influencing 
the qualitative behaviour of the solutions. Transient growth in the magnitude of the velocity 
components is found to be directly attributable to the compactness of the initial vorticity. 
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1 INTRODUCTION 


Vortex dynamics has been a significant subject of extensive investigation throughout the 
history of fluid mechanics. On the one hand, this interest has meant the actual establishment of 
the velocity compatible with a particular vortex with the results leading to numerous prototypical 
vortices that now carry names attributable to the person or the situation that created their 
existence. A Burgers vortex or a trailing vortex are but two of the numerous examples. Then, 
the question of the stability of any vortex has been asked. This is especially relevant if one 
wants to understand the consequences of a vortex under the influence of adverse or favorable 
strain. Finally, there is the topic of vortex breakdown, a special separate topic because of its 
consequences. Vorticity and vorticity interaction are the bases of turbulence. Heat and mass 
transfer rely to a large degree on vortex dynamics. And, there is the geophysical environment 
where vortices play major roles in both the ocean and the atmosphere. It is important, then, to 
know the behavior of vortices, particularly if they are subjected to external disturbances. 

The use of classical stability theory in examining a vortex flow has met with only limited 
success due to the fact that the mean flow is always a function of more than one independent 
variable. As a result, traditional normal mode analysis can only be used in an approximate 
manner. (This conclusion is always inevitable if the perturbation system cannot be reduced 
to one of ordinary differential equations.) Consequently, an alternative method for analysis 
must be sought in order to understand the dynamics. One such procedure is available and 
has already been applied to one or more rotating flows or vortices in a plane. Here, since the 
basic mean flow can be synthesized from linear functions of the appropriate spatial variables 
(constant shear), more general solutions of the perturbation quantities can be obtained than 
those derived from normal modes. In addition, there is the bonus that the perturbation field 
can be considered three-dimensional (cf. Craik and Criminale 1 for the basic method; Bayly 2 , 
Landman & Saffman 3 , Waleffe 4 , Dritschell 5 , or Craik and Allen 6 , among others, for specific 
applications of the method to planar vortical or rotating flows). Both approaches have found 
that such basic flows can be unstable and the instability depends upon various parameters and 
whether or not the perturbations are two or three-dimensional. Regardless of the conclusions 
from any of these treatments, it should be stressed that the work has relied on a basic vortical 
velocity that is two-dimensional. More realistic vortices have flow along the axis of rotation as 
well as motion in the plane. In other words, the mean flow is three-dimensional and investigation 
of the stability should consider this fact along with the possibility of the perturbations being 
three-dimensional. 

Any general treatment of the dynamics of perturbations in an arbitrary vortex is not possible 
analytically and is the principal reason why such work has been limited to or approximated in 
ways such as those cited above. Still, if one examines the family of vortices that is known and 
is willing to exchange one set of assumptions for others in order to obtain new insight, it is 
possible to investigate a three-dimensional vortex under three dimensional perturbations. This 
is the major point of this presentation. And, instead of seeking stability criteria per se, an 
initial- value problem is both posed and solved. The results are general, can be expressed in 
closed form, and both the early period transients as well as the asymptotic fate of disturbances 
can be determined explicitly. 


1 



2 BASIC PROBLEM 


One of the better references dealing with the derivation of vortical flow is that due to 
Donaldson and Sullivan 7 where a great deal of effort was given to deriving a general solution for 
a vortex family. It was found that the basic mean velocity can be written as 


u = (U r (r),V e (r),zW(r)) 


( 1 ) 


where the velocity components are given in the r, 9 , and ^-directions respectively. The coordi- 
nates (x,y) define the plane with r = ( x 2 + r/ 2 ) 1 / 2 , 9 = tan _1 (r//a;); z is perpendicular to the 
plane of rotation. This representation contains many of the more well-known examples, such as 
the Burgers’ vortex, Rott’s modified Burgers’ vortex, or the Batchelor trailing line vortex among 
others. In addition, (1) also contains viscous effects. The simplest fully three-dimensional vortex 
with this form of representation is the potential vortex under strain and is a limiting form of 
(!) or 

Ur = ~y, Vg = —, W = e. (2) 

2 r 

The coefficient To is proportional to the circulation, and the coefficient e, which can be time- 
dependent and even negative, is proportional to the strain rate. The perturbation problem 
associated with the vortex (2) can be solved explicitly as an initial- value problem. In effect, 
(2) is the inviscid limit of (1), and although being three-dimensional, it does not have a finite 
core. This lack in the basic flow does not carry over to the perturbation held. The pressure 
corresponding to (2) is 


P = - 


8 


(r 2 + 4 z 2 ) 


It 

2r 2 ' 


( 3 ) 


The vortex represented by (2) has zero vorticity. Thus, it is convenient to pose the pertur- 
bation problem in terms of the perturbation vorticity which is the leading term for the vorticity 
of the how itself. With the additional assumptions of an incompressible and inviscid huid, the 
perturbation equations for the vorticity are 


Du> r 

~dT 



= 0, 


and 


D^e 

Dt 


c 2r 0 

+ —Ue H o~ u r 

2 r 2 


= 0, 


Du z 

Dt 


— eu> z = 0 


where the differential operator is given by 


D d er d T 0 d d 

— = 1 1- ez — 

Dt dt 2 dr r 2 d6 dz' 


( 4 ) 

( 5 ) 

(6) 

( 7 ) 


The variables u r , ujg and uj z are the vorticity components in the r, 9 and z directions, respectively. 
In this differential operator, the spatial variables r and z appear linearly as coefficients of the 
partial derivatives with respect to r and z, allowing for a transformation of the independent 
variables to a moving coordinate frame given by 


£ 


= Mt)r, 


(8) 
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and 


£ = e io e ( t ) dt z = if) 2 {t)z, 


(9) 


r = /. (io) 

In this frame, advection due to the mean flow in these directions is removed from (4)-(6). At 
the initial time t = T = 0, the coordinates coincide as £ = r and ( = z, and there is no difficulty 
in applying initial conditions. Since the solution must be periodic in the 0 variable on physical 
grounds, the transformation of the dependent variables as / — ► fe m6 will be used throughout. 
On combining these two steps, the differential operator (7) reduces to 

_ _d_ inlp'ipl 

Dt dT e ’ 

and the governing set of partial differential equations now have coefficients that are functions of 
£ and T only. These equations are 



d(jJ r inTo^l e 

dT £ 2 2 ’ 

(12) 


dojg inT otpl e 2r O '0i 

QJ' + £2 tjJ ° ~ 2 Ue £ 2 Ur ’ 

(13) 

and 

duj z inToipl 

_ + UJ Z — €UJ Z . 

qji 1 c2 z z 

(14) 

The equations for uj r and uj z are uncoupled and can be solved directly, while the equation for 
ujg depends upon the solution for u r . On solving (12)-(14), it is found that 


LO r = U%£,C)4>l4>2E, 

(15) 


2 r i-T 

ue= ^(£,C)--^y o 1>l{T')dT'i4{U) M 2 E, 

(16) 

and 

u z =u° z {Z,01’lE 

(17) 

where = u r , ojg = 

: ujg, and uj z = uj z at T = 0. The exponential factor E is defined as 



E(£, T) = exp( — ^ £ ipl(T r ) dT'). 

(18) 


In all cases, \E\ = 1 and hence the factor E in no way contributes to the growth or decay of 
the amplitude of any quantity. It should be noted at this point, that if the straining parameter 
e is a positive constant, then the radial and tangential vorticity components will decay in time, 
while the maximum of the axial vorticity component will grow in time, with the radial location 
of this maximum approaching the core. This is certainly not surprising but it is the effects of 
this behaviour on the velocity held that will be of interest. 

The solutions (15)-(17) provide the means to describe the evolution of an arbitrary distur- 
bance to a potential vortex under strain. In order to do this, it is only necessary to prescribe 
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initial conditions in terms of the vorticity components subject to the constraint that the vortic- 
ity vector must be solenoidal. Thus, only two of the three components w)?(£, C), cj# (£,(), and 
4(£>0 are independent. At time T = 0, for example, two components are specified while the 
third is found from the relation 


10 (&*>?) , in ..o , 


94 

d( 


= o. 


(19) 


Once the vorticity has been determined, the perturbation velocities u = ( u r ,vg,w ) are derived 
by integrating the equations obtained from the definitions of the vorticity, namely 


and 


imjj-i dvg 

Ur = c W tp2 o a i 

£ d( 

(20) 

du r dw 

w ° = * 9( S( ■ 

(21) 

t/’i (d{£v g ) . \ 

= el mUr ) ' 

(22) 


The homogeneous solution to (20)-(22) represents the potential velocity held and is included in 
the integration in order to satisfy velocity boundary conditions. 


3 INITIAL PERTURBATIONS 

In this section, various initial conditions are examined in order to highlight the dynamics 
of the how with special attention given to the evolution of the velocity held. This is especially 
relevant in view of the fact that boundary conditions are in terms of the velocities rather than 
the vorticity. The actual solving of the partial differential equations for the velocity components 
can be avoided by using a Fourier transform in the ( variable without loss of generality. All 
dependent quantities can be transformed as 

/ OO 

mcv^dc. ( 23 ) 

-oo 

It is noted that only in the moving coordinate system will the Fourier transform provide the 
simplihcations that allow the analytical solution to progress. Equivalently, separable solutions in 
the moving coordinate can be sought, and this approach is taken here. Of course, initially two- 
dimensional disturbances can be considered and should the initial disturbance be independent 
of the ^-variable, for example, then for each n the governing equations are ordinary differential 
equations at the outset and the transform is not required. 


3.1 Impulsive Sources 

The first case considered is the two-dimensional disturbance with sources specified mathe- 
matically as Dirac delta functions in terms of the radial variable. This choice essentially provides 
a Green’s function for all arbitrary disturbances in this class. At time / = 0 and in the physical 
coordinates, the vorticity components are taken as 

4 = ~ r o), 4 = fte8'(r - ro), (24) 
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and, for consistency, it is required that 


(jj° r = —inSl, 


S(r - r 0 ) 
r 


(25) 


The full evolution of the vorticity is readily discerned from equations (15)-(17) while the axial 
velocity is found from (20) to be simply 


w = -Qff6(£ - r 0 )ip 2 E(£,T). 


(26) 


The radial velocity must satisfy 

^ < L^L + + (1 - n 2 )u r = -inQ, z £8(t - r 0 )ipiE(^,T) (27) 


which has the solution (for n > 0; analogous solutions for n < 0) 


_ i^r/’iT(r 0 ,T) J 

u r - 2 ] 

r 

u 4)"" 1 

o o 
£- 

V A 

uj> uj> 

(28) 

The corresponding tangential velocity is 




Q, z 'ip 1 E(r 0 ,T) j 
Ve= 2 { 

v r 0 7 

V r 0 7 

o o 

V A 

uj> uj> 

(29) 


To interpret the results, the velocities are transformed back to the physical coordinates r and t 
and the factor e m 6 is included. It might also be of interest to calculate another physical quantity, 
such as energy. For this example, both the velocities and the energy are determined. 

To highlight the differences of the strained vortex as opposed to one in the unstrained 
state, both e = 0 and e = eo 5 a constant, are considered. For the unstrained vortex, e = 0, 
ipi = ip 2 = 1? and in physical coordinates, E(r,t) = E(r,t)e m 6 = exp (in(9 — y§T)) with the 
velocity components becoming 


w = -Q, g 8(r - r 0 )E(r,t ), 


(30) 


and 


u r 


v e 


1^1 z E( tq , t ) 

2 


Q, z E(r 0 ,t) 

2 


' r \n — 1 
v r 0 ) 


' r 5—71 — 1 
v ro ) 


' r \n — 1 
v ro ) 


r 5—71—1 
r 0 ) 


r < r 0 
r > r 0 
r < r 0 
r > r 0 


(31) 


(32) 


It is noted that the velocity has a traveling wave solution in the angular direction. 

For constant straining, e = eo 5 'Pi = exp(yt), 'P 2 = exp( — Co t) and in moving coordinates 
E(£,t) = E(^i)e in6 = exp [in(6 — ~^i( eeot — 1)^). The resulting velocity components in phys- 
ical coordinates are 


u r 


W = —QgS( 

iQ z E(tq , t) 

2 


re 2 °‘- ro) exp( — eot)E(ripi,t), 
e 2 t ()^) n r < r 0 e 

< 

, r > r 0 e 


€Q_t 

2 


1 

€Q_t 


2 


(33) 

(34) 
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and 


v e = 


tt z E(r 0 ,t) 


t( r \n — 1 


— e 2 — 


_ TT~t( r \ — n — 1 


' r 0 ■ 


_ fo£ 
r < roe 2 


r > roe 2 


(35) 


For this case, there is a squeezing effect in which the overall pattern collapses into a core, the 
size of which is decaying exponentially in time. When £lg J 0 and Cl z = 0, the magnitude of 
the velocity decays and the traveling wave motion, described by the quantity E , decelerates. 
For Cl z 7 ^ 0 and fig = 0, there is an apparent exponential growth in the magnitude of the 
velocities u r and vg in the developing core. This is seen by noting that the velocities have their 

€Q t €0 

maximum values at £ = ro (or r = roe~~ ) and the maximum values grow as e~ . There is also 
a corresponding exponential increase in the velocity of the traveling wave. Thus, the two modes 
are seen to behave differently. It is unclear how these effects relate to the overall stability of the 
flow. 

To explore the dynamics from a different perspective, the disturbance energy density, given 

by 

1 f 00 

£n = 77 / (Kl 2 + \vg\ 2 + \w\ 2 )r dr (36) 

is conveniently calculated in the moving coordinates by 



■r | 2 + M 2 + M 2 V d£. 


(37) 


By considering only the case Q z J 0 and fig = 0, the energy is 


f - 


ny 0 

4 1 rz | 


(38) 


so that the perturbation energy for these two-dimensional disturbances remains constant during 
the entire evolution even though the velocity profiles are increasing in magnitude. This implies 
that the long time behaviour predicted by the Green’s function will most certainly be modi- 
fied when considering a spatially distributed source of vorticity. In the next subsection, this 
conjecture is shown to be true. 

Another set of initial conditions for which analytic solutions can be determined is taken as 
being axisymmetric (n = 0) with normal mode behaviour in the axial direction. Again, using 
Dirac delta functions to facilitate analytic solutions, the initial vorticity is specified as 


w r ° = n r S(r - r 0 )e-^ oZ 


(39) 


and 

co° = n e E(r-r o)e~^ z . (40) 

The value of V), which is not needed in order to determine the velocities, is determined by solving 
(19). The vorticity solutions of the initial value problem are given by equations (15)-(17) with 
the exponential factor E = 1. The angular velocity vg is found to be 

Vg = oCg(£ _ rQ y (41) 

To 
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To determine the radial and axial velocity components, it is noted that the solutions are separable 
in £ and £. Hence, by dropping the exponential factor e _ * 7o< “ in all quantities, the ordinary 
differential equations that the remaining velocities must satisfy are 


2 d 2 u r du 

£ J/-9 £ 


de 


dH 


1 + 


700 


2 f2 


U r = 


“*7002 & 


Si.-Si,^ 2 / T 1 fidT' 

£ Jo 


S(i - ro) (42) 


and 


The solutions of (42) and (43) are 


w = 


r/q d 
* 7 oV >2 d£ 


(£u r ) . 


u r 


-ilorp^l 

01 


ft, 


and 


= -''’ol/’f|7o 

0i 


ft, 


ft r ^ [ T J)\dT' 
r o Jo 


A'i(tl7oko)/i(J|7olO 

A'i(tl7o|£Ki(fel7oko) 


ft r ^ [ T J)\dT' 
r o Jo 


^r(tl7oko)/o(tl7olO 

A'o(Jl7o|£Ki(Sl7oko) 


(43) 


£ < j’o 
£ > j’o 


(44) 


£ < ro 
£ > ro 


(45) 


where ltd and /i are the modified Bessel functions. 

On examining the radial and axial velocities, it is found that in the absence of a straining 
vortex held, the terms proportional to ft$ are stationary in time whereas the magnitudes of the 
terms proportional to ft r grow linearly in T for all time. The addition of strain to the vortex 
held has a stabilizing effect in the following manner. The terms proportional to ft$ are always 
exponential in time and decay in magnitude. The terms proportional to ft r display algebraic 
growth for small time and the large time behaviour under constant strain, namely 


Uy* — y. 


2i7or 0 f 

eo r 0 


r 

ro 


e ~e 0 tro_ 

T 


_ 21 + 

r < roe 2 


_ 31 + 

r > roe 2 


(46) 


shows eventual exponential decay (the axial velocity component decays even faster). Even 
though standard stability based on large time behaviour would consider this how to be stable, the 
initial algebraic growth could be large enough so that nonlinear effects might become signihcant 
before the exponential decay dominates. It is stressed that the solutions are not restricted to 
constant strain. For example, if e(t) had the behaviour that e(0) 0 0 but e(t) — ► 0 for large 
time, the following scenario is possible: A period of algebraic growth is followed by a period 
exponential decline, and then the very large time behaviour returns to algebraic growth. Clearly, 
the overall effects depend upon the relative time scales. 

It is of quite some interest to compare the results of the axisymmetric initial disturbance 
with the two-dimensional initial disturbance. Although the vorticity components evolving from 
a two-dimensional initial disturbance have linear growth in time (which continues indefinitely for 
zero strain), the velocity components do not exhibit this behaviour. For this case, the term that 
produces the linear growth, namely dt' , only appears in the exponential factor that controls 

the traveling wave and not in the magnitude of the velocities. In the case of the axisymmetric 
disturbance, the linear growth term appears directly in the magnitude of the velocities. 
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3.2 Specific Initial Conditions 

In this section, special initial sources of vorticity are employed in order to understand how 
the vortex evolves in response to more general initial disturbances. Although complete analytical 
solutions for the velocities can be written in terms of the above Green’s functions, the resulting 
integrals are not easily analyzed. Therefore, numerical integration of the governing ordinary 
differential equations has been used instead. 

First, the two-dimensional disturbance (independent of the axial variable Q is considered 
and is taken to be 



= n z 

1 — ( r — tq ) 2 

: /I 6 ^ ’ 

V 47T (7 

(47) 

= 

—fie 

r — rn -0-r n ) 2 

P 4(7 

2<j\/4ircr ’ 

(48) 



~ r n -O-ro) 2 



infl { 

6 4(7 £ 4(7 

9 n 

(49) 



ry 47r<7 


These sources of vorticity are chosen such that, as a — ► 0, the above functions reduce to the Dirac 
delta function and its derivative centered at r$ and compare directly to the initial conditions 
(24), and (25). The velocities must satisfy 

~ r o -((- r o ) 2 

w = tt e ei \~ e " 1> 2 E(Z,T), (50) 

^V47T<7 

£ 2 (u r )tt + H(u r )i + (1 - n 2 )u r = U i<r 0) E(£,T), (51) 

V 47T (7 

and 

£ 2 ( v e)tt + 3£(^)£ + (1 - n 2 )v e = -^==— (Ve ^ E(^T))j (52) 

Equation (51) was solved numerically for various values of n and Co with a = .5, and fl z = 1.0. 
The results are shown in Figures 1 and 2 as functions of the real physical variables r and t. Figure 
1 gives the maximum magnitude of the velocity component u r as a function of time. Contrary 
to the prediction of the analytic solution of the Green’s function, the large time behaviour of 
this maximum is decay rather than exponential growth. The initial behaviour, on the other 
hand, does show transient growth that is dependent on the strain rate. Thus, for large strain 
rates, finite amplitude disturbances may cause the system to become nonlinear, even though 
the long-time solution is stable. In Figure 2, the velocity profiles for u r are shown for n = 2 
and e 0 = 1. Figures 2a-b are for an initial vorticity disturbance which is distributed in space 
(a = .5). This solution is equivalent to an integration of the Green’s function over ro, and 
produces a velocity profile that is decaying in time. Also to be noted is the oscillations as the 
disturbance is convected into the core by the positive strain. The origin of these fluctuations can 
be traced to the exponential factor and specifically it is a result of the spatial variation 

in the velocity of the traveling wave in the angular direction. It is interesting to note that the 
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0.0 1.5 3.0 0.0 0.4 0.8 1.2 


r (t = 2 ) r (t = 4) 

Figure 2: Plot of Re(u r ) (solid) and Im(u r ) (dashed) as a function of the physical coordinate 
r for spatially compact and spatially distributed initial sources of vorticity (fG = 1.0, fig = 0, 
e = 1.0, n = 2, To = 1.0 and ro = 2.0). Graphs for three values of a are shown: (a)-(b) a = .5, 
(c)-(d) a = .01, and (e)-(f) a = .001. 
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Figure 3: Plot of max\u r \ as a function of time for a spatially compact initial source of vorticity 
(fG = 1.0, fig = 0, o' = .001, To = 1.0 and ro = 2.0). Graphs for three values of n are shown: (a) 
n = 1, (b) n = 2, and (c) n = 4. The dashed curves are the Green’s function predictions. The 
symbols and the line through the symbols are the Green’s function prediction and the calculated 
response, respectively, for €q = —.2 
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Green’s function does not predict this behaviour directly since the exponential factor appears 
as E(ro,t). Therefore, the traveling wave velocity has no radial dependence. The velocity 
fluctuations and general decay in magnitude must be the result of the distributed nature of 
the initial source itself. To explore this phenomenon, Figures 2c-d and Figures 2e-f show the 
results of a numerical integration for a = .01 and a = .001, respectively. Both approximate the 
solution to a delta function source quite nicely for small times, as can be seen by examining 
the / = 2 graphs. The u r profile is linear near the origin and decays for r large as predicted by 
equation (34). The predicted response to the delta function initial condition continues until the 
cumulative effects of the exponential factor show themselves at / = 4. By comparison with the 
previous results, the rapid fluctuations are concentrated near the peak when a = .01 as opposed 
to the broader fluctuation when a = .5 and are not yet seen for a = .001. Of particular note 
is that the magnitude of the solution for a = .01 is decaying which is consistent with the prior 
results for the distributed source, and for the a = .001, the magnitude is still growing at / = 4 
which is consistent with the Green’s function results. The Green’s function predicts exponential 
growth for all time, but since there is a little bit of spatial variation in the initial profile, the 
a = .001 solution eventually goes the way of the a = .01 solution. This can be shown to be 
the case by examining the maximum of u r as a function of time for a = .001. The result of 
this calculation for the same values of n and Co as in Figure 1 are shown in Figure 3. For each 
case with Co > 0 there is an initial growth of the maximum, but then eventual decay. The 
decay is delayed for smaller values of Co which is consistent with the disturbance entering the 
core of the vortex as the source of the departure from the Green’s function predicted behaviour. 
Shown as dashed lines in Figure 3 is the growth predicted by the Green’s function solution. Of 
some interest is the behaviour for negative strain. Also shown in Figure 3 is the solution for 
Co = —-2 which decays and shows no initial growth. In this scenario, the initial disturbance is 
convected away from the core by the straining, and the Green’s function prediction shown as 
symbols is accurate. This investigation shows that the source of the fluctuations in velocity as 
the strain convects the disturbance toward the core are indeed a result of the spatial distribution 
of the initial source of vorticity and that this spatial distribution leads to the eventual decay of 
perturbations. 

The set of initial conditions corresponding to the axisymmetric modes (n=0) are taken as 


LO° = 4 l r 


roV 4 


-( r ~ r or 
-e 4cr e 


—27o z 


7T (7 


0 0 r - r 0 -(r-ro f ■ 

^ r . — e 4CT e 70 , 

LG V 47T G 

with lo ° obtained by requiring the vorticity to be solenoidal. The velocities satisfy 


Vg = - 


~ ( 7 r ° )2 / -hn C 
=e 4ct ipie , 
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d 2 u r ,du r 

? de dC 
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Tot/’ 2 
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7o r oV 4tt<7 
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Vh V4 vt(7 
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Figure 4: Plot of max\u r \ as a function of time for an axisymmetric (n = 0), spatially distributed 
initial source of vorticity with = 1.0, fig = 0, Tq = 1.0, a = .5, 70 = 1.0 and ro = 2.0. 
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with the last equation being integrated numerically. The maximum of u r as a function of time 
is shown in Figure 4 for various values of the strain rate. For this figure, fig = 0 and = 1. 
The above equation shows linear growth of the perturbation for the unstrained vortex e 0 = 0. 
For the strained vortex the same rate of linear growth is initially present and is independent of 
the choice of Co- However, as time progresses the exponential decay becomes evident with the 
decay occuring earlier for larger values of strain. The spatial profiles of the velocity are swept 
toward the center of the vortex by positive straining but, due to the absence of the traveling 
wave, the profiles remain smooth. For this case where the factor E = 1 and thus is not involved 
in an integral representation of the solution, the results of the Green’s function are predictive, 
with no further interpretation necessary. 


3.3 Three Dimensional Disturbances 

It has been noted that the quantity dt' plays a key role in the dynamics of two- 

dimensional disturbances (independent of the axial variable) and in the dynamics of axisym- 
metric disturbances. In the first instance this quantity controlled the speed of a traveling wave 
through the quantity Ee m 6 and noticeable effects are produced when considering a velocity pro- 
file that evolves from a spatially distributed initial source. In the second instance, the quantity 
Jq ipi dt' controlled the initial algebraic growth of the amplitude of the radial velocity compo- 
nent. For the fully three-dimensional disturbance, the relevant quantity should be FT m0_ *voC, 
which is a helical wave moving along the axial core. 

The form of the initial condition selected is given by (53) and (54) with the additional factor 
of e m6 . The vorticity is given by (15)-(17) after requiring the initial vorticity to be solenoidal 
in order to determine the uj z component. It is convenient to work with the equation for angular 
velocity given by 


r 2 dv e 

* de 


I V 0 I n 2 2 V-y f2i 

+ 3 ^^ + [f-ra -70712^ \ v e 


s 

>r 


1 d(tj 2 LO z ) 

i>! d £ 


, ' ^2 f2 
+ u r 

4>i 


(57) 


and the other velocity profiles are determined directly from the vorticity equations (20) and (22). 
The results of the numerical integration for the fully three-dimensional disturbances indeed show 
the artifacts from the previous special case analysis, but should be considered as a completely 
different class of solutions. In Figure 5, the amplitudes of the velocity are shown as a function 
of time for n = 2, ro = 4, 70 = 1, = 1.0 and ujg = 0.0. The dashed curves correspond to a 

delta-like initial condition with a = .01 and the solid curves show the evolution of a spatially 
distributed initial condition with a = .2. The radial velocity demonstrates exponential growth 
until the disturbance reaches the core of the vortex, and then there is rapid decline. As before, 
the spatially compact initial vorticity disturbance experiences the exponential growth for longer 
periods of time. The results for the other two velocity components are a little more ambiguous. 
For the tangential velocity, there appears to be a weaker exponential growth, especially for 
<7 = .01. For a = .2, there is very little growth of the tangential velocity. The axial velocity 
demonstrates basic exponential decay except for a short initial period where algebraic growth 
occurs. This is easily seen in Figure 5c where Co = 1-0. For Co = -2 there appears to be a balance 
between the algebraic growth and the slower exponential decay. 
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Figure 5: Plot of max\u r \, max\vg\ and max\w\ as a function of time for a spatially compact 
(a = . 01 , dashed) and spatially distributed (a = . 2 , solid) initial source of vorticity = 1 . 0 , 
fig = 0, n = 2, To = 1.0, 70 = 1.0, and r q = 4.0). Two values of the strain rate were chosen: 
(a-c) e 0 = 1-0 and (d-f) e 0 = 0 . 2 . 
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Figure 6: Plot of the real part (solid) and imaginary part (dashed) as a function of the moving 
coordinate £ for a spatially distributed initial source of vorticity (fi r = 1.0, fig = 0, a = .2, 
e = 1.0, n = 2, To = 1.0 and ro = 4.0). Graphs for three values of time t are shown: (a)-(b) 
t = 0, (c)-(d) t = 2, and (e)-(f) t = 4. 
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Figure 7: Plot of the real part (solid) and imaginary part (dashed) as a function of the moving 
coordinate £ for a spatially compact initial source of vorticity (fi r = 1.0, fig = 0, a = .01, 
e = 1.0, n = 2, To = 1.0 and ro = 4.0). Graphs for three values of time t are shown: (a)-(b) 
t = 2, (c)-(d) t = 4, and (e)-(f) t = 5. 
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In Figures 6 and 7, the angular and radial velocity profiles are shown as functions of the 
moving coordinate £ for Co = 1-0, r$ = 4.0, n = 2, 7 = 1.0 = 1.0, and fig = 0 which are the 

conditions for Figure 5a-c. The magnitudes of the velocities are increasing as time progresses. 
At around 1 = 4, the effects of the factor E start producing the fluctuations in the tangential 
velocity components with the fluctuations being more prevalent for the spatially distributed 
initial conditions than for the spatially compact initial conditions. For the case a = .2, the 
radial fluctuations in the radial velocity profile are not as visually evident, but have produced 
somewhat of a smoothing effect. For the case of a = .01 there are no discernible effects of 
the factor E on the radial velocity profile, even for 1 = 5. This is because the effects are still 
concentrated in a small interval around £ = 4 as can be seen by the figures for the tangential 
velocity. As time progresses, the solutions for a = .01 will show greater signs of the effects of 
the factor E and will decay in magnitude quite rapidly as can be determined from Figure 5. 

4 CONCLUSIONS 

The general problem of vortices subjected to perturbations is considered within a completely 
different framework than has heretofore been done. First, the nature of any vortex is reviewed 
with respect to the velocity field that might be required. Intrinsically, it is suggested that any 
vortex must be three-dimensional and such flows are usually under the influence of positive or 
negative strain. 

Second, the question of what might happen if perturbations are introduced into one of these 
flows has been a major challenge due to the fact that no direct mathematical analysis has been 
possible. The use of classical stability techniques that employ separable solutions is not possible 
in view of the underlying complications. Indeed, save for some numerical treatments related 
to special approximations, success has been only with those flows that are vortical in a plane. 
This work offers a means whereby these obstacles can be overcome and still include the full 
three-dimensionality of both the mean vortex and the perturbation field. 

If the strength of the vortex is strong enough and no boundaries are present, any effects 
due to viscosity are confined to the core. In effect, this means the core has a finite size and no 
component of the velocity is discontinuous. To a large degree of approximation, viscous effects 
can also be neglected in the analysis of the perturbation dynamics. Even on this basis, it is 
shown that the perturbation field does develop a core and no component of the velocity results 
in a discontinuity. And, the benefit of being able to obtain analytical closed form solutions for 
the initial- value, boundary- value problem is an asset. This has been possible for arbitrary initial 
values. Specifically, a potential three-dimensional vortex is subjected to perturbations at time 
t =0 and the full three-dimensional problem is solved for all values of time. 

The strategy used for solution directly involves the perturbation vorticity. The linearized 
equations governing the vorticity are solved subject to initial distributions specified as functions 
of the space variables. The kinematic requirement that the vorticity be solenoidal is insured. 
Then, these results are used to solve for the velocity and meet proper boundary conditions. For 
localized vorticity (specified as Dirac delta functions), there is transient growth in the magnitudes 
of the velocity components. It is also found that the larger the strain rate, the larger the transient 
growth and this goes on until the disturbance enters the core of the vortex. For distributions 
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of the perturbation vorticity that are more evenly distributed in space, there is little, if any, 
transient growth for the velocity. Asymptotically, no disturbance survives regardless of whether 
or not it is localized. The special cases of a two-dimensional or an axisymmetric vortex are 
included. Even though diffusion is absent in this work, it appears that, once vorticity is diffuse 
- as can be inferred from the more evenly distributed initial input - then even severe strain 
cannot enhance the flow. 
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